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Summary. We study the competition between several advantageous mutants in 
an asexual population (clonal interference) as a function of the time between the 
appearance of the mutants At, their selective advantages, and the rate of deleterious 
mutations. We find that the overall probability of fixation (the probability that at 
least one of the mutants becomes the ancestor of the entire population) does not 
depend on the time interval between the appearance of these mutants, and equals 
the probability that a genotype bearing all of these mutations reaches fixation. 
This result holds also in the presence of deleterious mutations, and for an arbitrary 
number of competing mutants. We also show that if mutations interfere, an increase 
in the mean number of fixation events is associated with a decrease in the expected 
fitness gain of the population. 

1 Introduction 

Evolution, according to Dobzhansky (1973), is the unifying concept that pulls 
together all the different strands of biology. Indeed, while evolution provides 
the framework to understand the otherwise bewildering panoply of adapted 
forms, it also allows us to understand the patterns that mutation and se- 
lection leave in the molecules of life, namely DNA and proteins. One of the 
central problems of the branch of biology that is devoted to the systemat- 
ics of the tree of all living things, molecular evolution, concerns the rate of 
adaptation of individual organisms and species. A precise understanding of 
molecular phylogeny requires accurate models of evolution and adaptation. 
In this contribution, we model the adaptation of chromosomes that do not 
undergo recombination, or more generally, the rate of adaptation of asexual 
organisms. 

The main observation that we study here is that the rate of adaptation of 
asexual organisms, or of regions of low recombination in the genomes of sexual 
organisms, does not steadily increase with increasing rate of advantageous 
mutations. When two or more advantageous mutations appear in different 
organisms at approximately the same time, then only one of them can go 
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to fixation, i.e., become shared by all members of the population, while the 
others will be lost. This loss of potentially beneficial mutations limits the rate 
of adaptation to the rate at which individual mutants can go to fixation. By 
contrast, in sexual organisms, several mutations can recombine and thus go to 
fixation together. This interference of advantageous mutations has long been 
recognized as a potential disadvantage of asexual populations (Fisher, 1930; 
Muller, 1964; Hill and Robertson, 1966). Recently, several groups have worked 
on an exact quantification of the interference effect, both in theoretical studies 
(Barton, 1995; Gerrish and Lcnski, 1998; Orr, 2000; Gcrrish, 2001; McVean 
and Charlesworth, 2000) and in experimental studies with bacteria (de Visser 
et al., 1999; Rozen et al., 2002; Shaver et al., 2002) and viruses (Miralles 
et al., 1999; Cuevas et al., 2002). A good quantitative understanding of the 
interference effect is necessary in order to assess the influence that interference 
has on the patterns of molecular evolution and variation in large populations. 

There are two separate dynamics that both contribute to the overall in- 
terference effect: First, if two advantageous mutants are both present in suffi- 
ciently high concentrations, such that loss to drift can be neglected, then they 
compete deterministically, and the mutant with the higher selective advantage 
will replace the other one. Second, if at least one advantageous mutant is still 
very rare, then we have to consider the influence of other mutants' presence on 
the chance that this mutant is lost to drift. To date, there is no single theory 
that takes into account both dynamics to their full extent. Gerrish and Lenski 
(1998) calculated the speed of adaptation as a function of population size and 
beneficial mutation rate, under the assumption that clonal interference can be 
neglected during the initial phase of drift. They assumed that the probability 
that a mutant is not lost to drift corresponds to one minus the standard prob- 
ability of fixation (as calculated for example by Fisher 1922, Haldane 1927, 
Kimura 1962) of that mutant. Orr (2000) modified the calculations of Ger- 
rish and Lcnski to include beneficial mutations that arise in genomes bearing 
one or more deleterious mutations, but also did not address the effect of in- 
terference of other beneficial mutations during the initial phase of drift. In 
general, the influence of deleterious mutations on the probability of fixation 
has been studied extensively (Manning and Thompson, 1984; Charlesworth, 
1994; Peck, 1994; Johnson and Barton, 2002), but there are very few studies 
that consider the effect of interfering beneficial mutations. The problem of 
interfering advantageous mutations is that there exists no simple theoretical 
framework with which their influence on drift can be described accurately. 
Barton (1995) derived an approximation that allowed him to calculate the 
probability of fixation of a mutant in a population that is undergoing a selec- 
tive sweep. However, his approximation is valid only for very small selective 
advantages. McVean and Charlesworth (2000) studied a similar situation with 
numerical simulations, and applied their results to codon bias and levels of 
polymorphism in molecular evolution. 

Here, we study the mutual interference of two advantageous mutants that 
are each initially present in only a single organism. We derive a phenomcno- 
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logical description of the mutants' fixation probabilities in the interference 
regime. This phcnomcnological description, which is in essence an interpola- 
tion between limiting cases that can be described with standard theory, agrees 
very well with numerical simulations. We also find that the expected fitness 
increase of the population is maximized if the mutant with higher selective ad- 
vantage arrives earlier, even though this order of appearance leads, on average, 
to fewer mutants that go to fixation. 

2 Model 

We consider three distinct sequence types: wild type, advantageous mutant 
1, and advantageous mutant 2, with fitness values 1, 1 + s\, and 1 + ,S2, 
respectively. Thus, si denotes the selective advantage of mutants of type 1, and 
S2 is the selective advantage of mutants of type 2. For simplicity, we assume 
that S2 > s\. Initially, the population is homogeneous and consists of the wild 
type only. Then, one randomly chosen wild type individual is replaced by an 
advantageous mutant of either type, and some At generations later, another 
randomly chosen individual of the population (cither wild type or mutant) 
is replaced by a mutant of the other type. Throughout this paper, we will 
understand At to be the difference in generations between the appearance of 
mutant 2 and mutant 1, so that negative values of At indicate that mutant 2 
appeared before mutant 1. 

The population is finite of size N, and replication takes place in discrete 
generations, according to the Wright-Fisher model: All individuals in genera- 
tion t are direct descendants of the individuals of the previous generation; the 
probability that an individual is the offspring of a particular parent is pro- 
portional to the parent's fitness. We introduce deleterious mutations into the 
offspring organisms with probability u. If a mutant of type 1 or 2 is hit by a 
mutation, then its fitness is set to 1 (that is, it reverts to the lower fitness wild 
type). We do not consider deleterious mutations in the wild type genotype. 

We consider a genotype to be fixed if it has become the most-recent com- 
mon ancestor of the whole population, regardless of whether some individuals 
in the population have a different genotype. This definition of fixation has 
been used recently to study fixation of beneficial mutations in a heteroge- 
neous genetic background (Barton, 1995; Johnson and Barton, 2002), and 
also to investigate the process of fixation in a viral quasispecies (Wilke, 2003). 

3 Theoretical analysis 

3.1 Probability of ultimate fixation 

In the Appendix, we derive equations for the probability of fixation P(s,u) and 
time to fixation T(s, u) of an individual mutant with selective advantage s and 
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mutation rate it. For two mutants, we are mostly interested in the probability 
of ultimate fixation, that is, the probability that a mutant reaches fixation 
and is not subsequently replaced by the other mutant. In the following, we 
will denote the probability that mutant i reaches ultimate fixation under the 
condition that the two mutants are introduced At generations apart as TTi(At). 
We can calculate ni(At) for the two limiting cases At — ► — oo and At — > oo, 
and present a phenomenological description for intermediate At. 

For At — > — oo, that is, when mutant 2 arises much earlier than mutant 1, 
we have 

7ri(-oo) = [l-P(* 2 ,«)]P(«i,u), (1) 

7T 2 (-C0) = P(s 2 , U ). (2) 

Since s 2 > Si, mutant 1 can reach fixation only when mutant 2 has not reached 
fixation. For At — > oo, we have to consider the possibility that mutant 1 
reaches fixation first, but is later replaced by mutant 2. Therefore, we find 

*i(oo) = pm [i - ^( (1+ 1 3 +; 2 _ M) - 1,«)] , (3) 

tt 2 (oo) = P(s u u)p( (1+ 1 g ^ ( S 1 2 _ u) ~ + [1 - P(s 1 ,u)]P(s 2 ,u) , (4) 

where P( (i + l^yi- u ) ~ is the probability that mutant 2 goes to fixation 
after mutant 1 has already reached fixation. Since P( — l, u ) is 

always smaller than or equal to P{s 2l u), we have 7Ti(— oo) < 7Ti(oo) and 
~k 2 {— oo) > 7T2(oo). In other words, both mutants have a higher probability of 
fixation when they are introduced first than when they are introduced second. 

For intermediate values of At, there is no theory that enables us to de- 
rive a simple expression for TTi(At) (but see Barton 1995). We cannot use 
branching process theory, because it assumes that the presence of the invading 
mutants does not influence the mean fitness (which they do for intermediate 
At). Also, diffusion theory becomes unwieldy when there are more than 2 
different sequence types. Nevertheless, we have found that we can develop a 
phenomenological description of the competition of two mutants as follows. 
We know that 7Tj(i) must reach the two limiting values 7Ti(oo) and 7Tj(— oo) 
for sufficiently large positive or negative At. Moreover, as long as At < 0, 
we do not expect iTi(At) to be very different from 7Tj(— 00), because mutant 

1 has only a realistic chance of proliferating and going to fixation if mutant 

2 is not present in the population. As long as At < 0, mutant 2 will either 
go to fixation relatively unscathed from the later appearance of mutant 1, 
or it will be lost to drift, in which case mutant 1 will have its turn. Like- 
wise, for At larger than the time to fixation of mutant 1, T\, we expect that 
TTi(At) — 7Tj(oo), because T\ generations after the introduction of mutant 1, 
mutant 2 will either find a population in which mutant 1 has already gone 
to fixation, or one in which it has been lost to drift. For < At < Ti, we 
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expect that ir 1 (At) smoothly increases from 7Ti(— oo) to 7r 1 (oo), while 7r 2 (Z\t) 
smoothly decreases from 7r 2 (— oo) to 7r 2 (oo). These considerations suggest a 
sigmoidal form for iTi(At). We use a logistic growth model to describe the 
smooth transition from 7Tj(— oo) to 7Tj(oo) within the range < At < T\\ 

. . . . . 7T,(oo) — 7T,-( — Oo) 

MAt) = M-oo) + ^^trj) > ( 5 ) 

where 71 is the fitness advantage of mutant 1 at finite mutation rate (see 
the definition following Eq. (8) below). This expression appears to be an 
acceptable description of the exact time dependence (see Numerical Results) 
without free parameters. 



3.2 Overall probability of fixation 

Besides the individual probabilities TTi(At) and 7r 2 (Z\i), their sum 7r = 
ni(At) + TT2(At) is also of interest. This sum is the overall probability that at 
least one mutant goes to fixation. When we sum Eqs. (1) and (2), or Eqs. (3) 
and (4), we find that the overall fixation probability does not depend on the 
order in which the mutants are introduced, and has the form 

7T = P(si, u) + P{s 2l u) - P( Sl ,u)P{s 2 ,u) . (6) 

Furthermore, if Eq. (5) is indeed a good description of iri(At) for arbitrary 
At, then it should not depend on At at all, because all time dependencies 
cancel when we use Eq. (5) to calculate -K\{At) + 7r 2 (Z\t). 

The invariance of the overall fixation probability n under the order by 
which the mutants appear is a general property. It holds also when more than 
2 mutants arise: Assume that n advantageous mutants arise, with selective 
advantages s\,...,s n compared to the wild type. Further assume that the 
time intervals between the appearances of the mutants are large. Then, the 
probability that none of the mutants make it to fixation is n"=i[l ~ P{ s ii u )]i 
regardless of the order of their appearance. The probability that at least one 
mutant goes to fixation is therefore 

n 

7r = l-JJ[l-P(s i ,«)] ) (7) 

i=l 

which reduces to Eq. (6) in the case of n = 2. Using an extension of Kimura's 
well-known result for P(si,u) derived in the Appendix, Eq. (A. 4), we can 
simplify this expression even further. We find in the limit N — » 00: 

7T = 1 - cxp [-2 ( 7l + 72 + • • • + 7 „)] , (8) 

with 7i = (1 + Si)(l — u) — 1 for u < Si/(1 + Si), and 7 j = otherwise. 
Moreover, in the absence of deleterious mutations, for u = 0, tt becomes 
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tt = 1 - cxp [-2 (si + s 2 + • •• + «„)] . (9) 

In this limit, the overall probability of fixation is the same as the probability 
of fixation of a single mutant with selective advantage s equal to the sum of 
the selective advantages of all invading mutants, s — Sj. Such a hypothet- 
ical mutant is extremely unlikely in clonal populations, because all beneficial 
mutations would have to hit the lineage sequentially, but could occur when 
beneficial mutations are shared via recombination. In the limit that all ji 
vanish (that is, for very large mutation rates or for vanishing si), Eq. (7) 
becomes 



3.3 Expected fitness increase and expected number of fixed 
mutants 

Depending on which mutant goes to fixation, the average fitness of the final 
population is either 1, I+71, or I + 72, with 7 as defined following Eq. (8). The 
expected fitness increase (j(At)) after the introduction of the two mutants is 
therefore 

(-f(At)) = 7i7Tl(^t) + 727T 2 (Zit) . (11) 

Since we know that ni(At) + 7r 2 (Z\t) is constant, and 72 > 71 as a direct 
consequence of our assumption s 2 > Si, it follows that (7(— 00)) > (7(00)). If 
wc write 7Ti(oo) = 7Ti(— 00) + Aw, and 7r 2 (— 00) = 7r 2 (oo) + An, then we find 

(7(-oo))-<7(oo)) = (72-7i)^r>0. (12) 

Thus, the expected fitness increase is larger if we introduce the mutant with 
the higher selective advantage first. 

Let 7ifi x denote the number of fixed mutants, that is, we have rifi x = 
if none of the mutants reach fixation, rifi x = 1 if exactly one of the mutants 
reaches fixation, and rifi x = 2 if first mutant 1 reaches fixation and is later 
replaced by mutant 2. For At — > — 00 and N large, na x can never be larger than 
one, because the probability that mutant 1 goes to fixation in the background 
of mutant 2 is zero. We find in this limit for the expected value of n^: 

(nfix (-00)) = 7ri(-oo) +7r 2 (-oo) = it. (13) 

In the limit At — > 00, on the other hand, we find 

(n fix (oo)) = tti(-oo) + 7r 2 (-oo) + P(s!,u)p( + \~^^_ M ) - 

= - + P(s,u)P( {1+ 1 s ^_ u) -l,u). (14) 

Clearly, (nfi x (oo)) > (rifi x (— 00)). This means that if the mutant with the 
smaller selective advantage appears before the mutant with the larger selective 
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advantage, then the expected number of mutants that go to fixation is larger 
than if the mutant with the larger selective advantage appears first. However, 
at the same time the expected increase in average fitness is smaller, as we saw 
in the previous paragraph. 

4 Numerical Simulation 

In order to measure fixation probabilities, we carried out 100,000 replicates 
of the simulation for each set of parameters, and recorded the final outcome 
(all individuals unmarked, or marked as descendants of either advantageous 
mutant). We studied a population of size N = 1000 and mutation rates u = 
0.0, 0.01, 0.02, 0.03, 0.05, 0.07, 0.1, 0.15, 0.20, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8. The 
selective advantages were si = 0.1 and s 2 = 0.2, s\ = 0.1 and s 2 = 0.5, 
si = 0.05 and s 2 = 0.2. We also studied a population of size N = 10, 000 
for a subset of these parameters, in order to make sure that our results were 
robust against a change in population size. Because of the sizable amount of 
CPU time needed to carry out 100,000 replicates for N — 10,000, we could 
however not study this case exhaustively. 

In order to keep track of the evolutionary history of each mutant, we 
marked the initial mutants 1 and 2 with two distinct inheritable neutral mark- 
ers. In that way, we could distinguish wild type sequences that were descen- 
dants from mutants 1 or 2 from the wild type sequences that were originally 
present. We continued all simulations until all individuals in the population 
were either unmarked, marked as descendants of mutant 1, or marked as de- 
scendants of mutant 2. 

Figure 1 shows the probability of ultimate fixation of mutants 1 and 2, m 
and 7T2, and the overall probability of fixation n — m + 7r 2 , for u = 0. We see 
that the probability n\ is constant for At < 0, and starts to increase as soon 
as At turns positive. Eventually, m levels off again. We observe the opposite 
behavior for the probability 7r 2 . For At < 0, 7r 2 is constant, but decreases 
rapidly in the same range of At in which m increases. Finally, 7r 2 levels off as 
well. 

The overall probability of fixation tt is approximately constant for all values 
of At. The solid lines in Fig. la and b correspond to Eq. (5), and the solid line 
in Fig. lc is m(— oo)+7r 2 (— oo). We find that our phenomenological description 
Eq. (5) performs very well for intermediate At. 

In Fig. 2 we plot the same quantities as those shown in Figure 1, but now 
with a positive mutation rate u — 0.05. We observe that the two probabilities 
7Ti and 7r 2 have smaller values than in the absence of mutations. The phe- 
nomenological description still works well, and appears to correctly take into 
account the effects of mutation. The overall fixation probability 7r is again 
independent of At. In Figure 3 we show the time to fixation of mutant 1, T\, 
and the time to fixation of mutant 2, T 2 , as functions of At. The parameter 
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Fig. 1. Probability of fixation as a function of time interval At (N = 1000, si = 0.1, 
S2 — 0.5, u — 0.0). Solid lines represent the prediction according to the logistic 
growth model Eq. (5). a: Probability of fixation of genotype 1 (less beneficial muta- 
tion), 7ri. b: Probability of fixation of genotype 2 (more beneficial mutation), 7T2. c: 
Overall probability of fixation, n. 



values are the same as those of Figure 1. In Figure 3a, we see that T\ is ap- 
proximately constant for all values of At except for the range < At < 200. 
For negative At the fixation of mutant 1 occurs only when mutant 2 has been 
eliminated. Therefore, T\ corresponds to the result for the fixation time of a 
mutant with selective advantage s\ in a homogeneous population with wild- 
type individuals only. The same is true for large positive At because there 
mutant 1 has enough time to reach fixation without the interference of the 
second mutant. For small positive At, we observe a decrease of T\. The de- 
crease occurs for those At for which the two beneficial mutants coexist for 
several generations in the population. Fixation of mutant 1 in this regime 
occurs only when mutant 1 reaches fixation so quickly that mutant 2 has not 
had time to build up momentum. Otherwise, most likely mutant 1 will be 
displaced by mutant 2 before reaching fixation. For mutant 2, the time to 
fixation T 2 is shorter for negative At than for positive At. For positive At, in 
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Fig. 2. Probability of fixation as a function of time interval At (N = 1000, si = 0.1, 
S2 — 0.5, u — 0.05). Solid lines represent the prediction according to the logistic 
growth model Eq. (5). a: Probability of fixation of genotype 1, 7Ti. b: Probability of 
fixation of genotype 2, TT2. c: Overall probability of fixation, it. 



a fraction of cases mutant 2 has to go to fixation in a background of mutant 
1, rather than in a background of wild type. In the background of mutant 1, 
the selective advantage of mutant 2 is smaller than in a wild-type population, 
which explains the increased time to fixation. Interestingly, before T 2 starts 
to rise for increasing positive At, it quickly spikes at At « 30. This spike 
has the following explanation: If mutant 1 is introduced right before mutant 
2 has reached fixation, then the time to fixation of mutant 2 is increased by 
the additional time it takes for mutant 1 to disappear again. This additional 
time will typically be one or two generations, which is in agreement with the 
height of the spike. (7\ on the top of the spike is actually approximately half 
a generation larger than before the spike.) The limiting values for m and 7r 2 
when At — > — oo and At — > — oo were estimated by means of branching pro- 
cess theory which provides very good accuracy, especially for low mutation 
values. 
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Fig. 3. Time to fixation as a function of the time interval between beneficial 
mutants At (TV = 1000, si = 0.1, S2 = 0.5, u = 0.0). a: Time to fixation of genotype 
1, Ti. b: Time to fixation of genotype 2, T2. 
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Fig. 4. Time to fixation as a function of the mutation probability u (si =0.1 
and s 2 = 0.5, At = 100, N = 1000). Solid line: the time Ti to fixation for the less 
beneficial mutant, dashed line: time to fixation of more beneficial mutant T2. 



In Figure 4 we show the fixation times Ti and T 2 as a function of mutation 
rate obtained from simulations, and compared to the prediction from diffusion 
theory, Eq. (A. 5). As we increase the mutation rate u, the population moves 
from a strong selection regime, characterized by a short time to fixation, to a 
neutral regime, where T w 2N. We also observe that the transition between 
these two regimes occurs at different mutation rates for the two genotypes: T\ 
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Fig. 5. Overall probability of fixation as a function of mutation rate u. Circles 
correspond to the situation of two invading mutants (s\ = 0.1, S2 = 0.5, At = 100, 
N = 1000), and squares correspond to the situation of three invading mutants 
(si = 0.8, S2 = 0.9, S3 = 1.0, At = 20). The solid line is the prediction according to 
diffusion theory, while the long-dashed line corresponds to the branching process. 
The short dashes indicate the fixation probability levels according to the neutral 
theory Eq. (10), for two or three mutants, respectively. 



shows an abrupt transition around the critical value u = U\ c « 0.09, whereas 
T 2 reaches the neutral regime around u = u 2c ~ 0.30. The mutation rate at 
which this transition occurs is known as the error threshold (Eigen, 1971). 

In Figure 5 we can see the overall probability of fixation as a function 
of mutation rate for the case of two as well as three interfering mutants. 
Simulation results are again compared to the diffusion theory result Eq. (A. 4), 
but also to a prediction from branching process theory, Eq. (A. 14). For these 
simulations, we used s\ — 0.1, s 2 = 0.5 and population size N = 1000. 
We also chose a time interval At = 100, which is smaller than the time 
required for fixation of mutant 1, to ensure that the dynamics takes place 
in the clonal interference regime. Above the error threshold, the probability 
of fixation is known to be 1/N for a single mutation (2/N and 3/N for two 
or three mutants, respectively). Because the branching process description 
assumes infinite population size, it predicts zero probability of fixation above 
the error threshold. Diffusion theory, on the other hand, describes this regime 
adequately, while being less accurate at small mutation rates. The abrupt 
change between the ordered and disordered regime predicted by both theories 
is not present in the numerical data because the simulated system is finite. 
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Fig. 6. a: Expected increase in fitness due to introduction of mutants as a function 
of time between their introduction, b: Expected number of fixed mutants (panel b) 
as a function of At (si = 0.1, S2 = 0.5, u = 0.0, N = 1000, for a and batt). Solid 
line in panel (a) is expected result according to fixation probabilities given by Eq. 
(5). 



Finally, we tested the prediction that the mean number of fixations in the 
population is larger when the mutant with the smaller selective advantage 
is introduced first, even though the mean fitness increase is smaller. Figure 
6a shows the results for the expected increase in fitness as a function of the 
time interval At as defined in Eq. (11), as well as the solution of Eq. (11) 
using a sigmoidal ansatz for the probabilities tti and 7r 2 , as defined in Eq. (5). 
As expected, the mean fitness gain decreases as At turns positive, while the 
expected number of fixed mutants, shown in Figure 6b, increases. 



5 Discussion 



The rate of adaptation of a non-recombining chromosome or asexual popula- 
tion is a non-trivial quantity that depends on details of the fitness distribution 
of mutations. Beneficial mutations often take a long time to dominate in a 
population, and compete against each other in the background of deleterious 
mutations. The probability of fixation of any one mutation depends on the 
time of their introduction, stochastic drift, and the prevalent mutation rate. 
A good understanding of the rate of adaptation of asexuals is important in 
molecular phylogeny, because it influences the speed of the assumed molec- 
ular clock (Kuhner and Felsenstcin, 1994). Since successive fixation events 
determine the branchings in the phylogenetic tree, a stochastic analysis of the 
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fixation probability of competing mutations can provide insight into models 
of evolution used for tree reconstruction methods (Page and Holmes, 1998). 

We studied the probability of fixation of beneficial mutations in the pres- 
ence of other beneficial mutations that were introduced either earlier or later, 
and in the presence of other deleterious mutations, that is, we studied fixation 
in a non-equilibrium background. This problem has been addressed previously 
by Barton (1995) using a deterministic approach, which was unfortunately 
limited to small selective advantages. Similarly, Johnson and Barton (2002) 
studied the probability of fixation in a changing background, but they con- 
fined themselves mostly to the case where deleterious mutations accumulate 
in the background after a selective sweep. 

We found that the probability of fixation is easily understood from an 
analytic point of view as long as the time interval between the introduction 
of mutants is larger than the fixation time of the least beneficial mutation. 
In the interference regime, where the population consists of clones of the wild 
type as well as both beneficial mutants, a phenomenological approach based 
on logistic growth is successful at describing the competition. In general, we 
find that the probability of fixation of the most beneficial mutant is reduced 
when it competes against a less beneficial mutant, and there is a sizable prob- 
ability that the less beneficial mutant will in fact survive (see Fig. 1). Yet, the 
probability that either one of the mutants survives does not depend on the 
timing of their introduction. 

A number of authors noted the reduction of fixation probabilities caused by 
segregating deleterious mutations (Peck, 1994; Manning and Thompson, 1984; 
Charlesworth, 1994; Orr, 2000; Johnson and Barton, 2002), and we observe 
the same qualitative behavior here. We find that a background of deleterious 
mutations suppresses the fixation probability, and this suppression is more 
pronounced for the less beneficial mutant. In the model presented here, there 
is a simple relationship between beneficial and deleterious mutations: Each 
mutation occurring on a clone created by a beneficial mutation is reduced to 
the wild type fitness. This simplified model corresponds to a single-peak land- 
scape, exhaustively studied in the context of quasispecies theory (Tarazona, 
1992; Swetina and Schuster, 1982; Galluccio, 1997; Campos and Fontanari, 
1998). Consequently, we expect a more pronounced decrease of the probabil- 
ity of fixation for small s due to the existence of an error-threshold, i.e., a 
point at which information is lost from the sequence due to a critically high 
mutation rate, see Eigen (1971). In the stochastic regime of evolution, which 
occurs for mutation rates above the error threshold, individuals replicate ran- 
domly, and each member of the population is expected to be the most recent 
ancestor of the population with the same probability n — 1/N while the time 
to fixation is T = 2N (Kimura, 1968; Kingman, 1982; Donnelly and Tavare, 
1995). This prediction is confirmed in Figure 4, where we see that for large u 
both Ti and T 2 are approximately equal to 2N. Accordingly, in the interval 
uic < u < U2c, where Ui C denotes the error-threshold for genotype of type i, 
we can describe the dynamics by ignoring the less beneficial mutant. 
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The explicit modeling of stochastic interference between beneficial mu- 
tation revealed an interesting observation concerning the mean number of 
fixation events and the expected fitness change. Naively one might surmise 
that, as the number of fixation events increases, so does the expected increase 
in fitness. Instead, we found the opposite dynamic: Fig. 6 clearly shows that if 
the mutant with the smaller selective advantage appears first, then we expect 
more fixation events because the second mutant then has a good chance to 
survive. Yet, this is not the best scenario if we are interested in maximizing 
fitness: In this case, it is better if the most beneficial mutant can establish 
itself first, even though this implies that there is little chance for another 
fixation event. 

Finally, we observed that the probability for any beneficial mutant to go 
to fixation is equal to the probability that a genotype bearing all mutations 
will become the ancestor of the entire population. Moreover, this observation 
holds for an arbitrary number of competing mutations, even in the presence 
of deleterious mutations, independently of the order of appearance of the ben- 
eficial mutations. Naturally, while such a probability is very high if the sum 
of all benefits is large, the probability of a single sequence appearing with this 
combination of mutations is exponentially small in non-recombining popula- 
tions. This result was confirmed using an extension of Kimura's result (1962) 
to finite mutation rates [Eq. (A. 4)] and also by means of a branching pro- 
cess formulation. Both approaches yield similar results and are in excellent 
agreement with our numerical simulations. We find that the branching pro- 
cess approach is more accurate at low mutation rates, away from the error 
threshold. As it is formulated in the N — > oo limit, it fails to predict the 
asymptotic result 7r s=s n/N which is reached in the limit that all 7$ vanish. 

The present analysis concerns the rate of adaptation as beneficial muta- 
tions compete against each other in the presence of deleterious mutations, in 
a non-equilibrium framework. This is a step towards the goal of characterizing 
the rate of adaptation of whole populations in arbitrary circumstances, and in 
the presence of recombination. We expect that these extensions are necessary 
before accurate predictions of optimal mutation rates in biological evolution 
can be made. 
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A Probability of fixation and time to fixation for an 
individual mutant 

We can calculate the probability of fixation of an individual mutant either 
from branching process theory (Fisher, 1922; Haldane, 1927; Fisher, 1930; 
Barton, 1995; Johnson and Barton, 2002; Wilke, 2003) or from diffusion the- 
ory (Kimura, 1962). Branching process theory is applicable for very large 
population sizes, and arbitrary (but positive) selective advantages s. Diffusion 
theory is applicable to moderately large to large population sizes and small 
or even vanishing selective advantages s. Moreover, diffusion theory also gives 
an expression for the expected time to fixation. 

A.l Diffusion theory 

If the advantageous mutant suffers additional mutations while it goes to fix- 
ation, then we typically have to use multi-dimensional diffusion equations, 
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which can be very unwieldy. However, in the simple case treating the advan- 
tageous mutant and the wild type only and where the only effect of deleteri- 
ous mutations is to revert the advantageous mutant back to wild type, one- 
dimensional diffusion theory, as developed by Kimura (1962); Kimura and 
Ohta (1969), is still applicable, and gives satisfying results (van Nimwegen 
et al., 1999; Wilke et al., 2001). 

The main quantities in diffusion theory are the mean M$ x and the variance 
Vg x in the rate of change per generation of the concentration x of the invading 
mutant. The derivation of these quantities is straightforward (Ewens, 1979), 
and we find 

M Sx = jx(l - x) , (A.l) 
V Sx = x(l - x)/N (A.2) 

to first order in 1/N and 7. Here, 7 is the relative difference in mean fitness 
between a wild-type population and a population with fixed advantageous 
mutant. The only difference between these expressions and those of Kimura 
(1962) is that 7 replaces the selective advantage s in M§ x . Without mutations 
(u = 0), 7 is equal to s, and we recover the standard results. For positive u, 
7 has the form 

f(l + *)(l-u)-l for «<«/(! + «), / a o\ 

7 \ fort* >«/(! + *)■ 

The mutation rate u c at which 7 reaches 0, u c = s/(l + s), corresponds to 
the error threshold of quasispecies theory (Swetina and Schuster, 1982; Eigen 
et al., 1988, 1989). For mutation rates larger than u c , the invading mutant does 
not have an advantage over the wild type, and the fixation process corresponds 
to that of a neutral mutant. 

The probability of fixation follows now from Kimura (1962) as 

P M = ^f^W ( A - 4 ) 

Likewise, the expected time to fixation follows from Kimura and Ohta (1969) 

as 

T(s,u) = J 1 + 1 ~J {s \ u) J 2 , (A.5) 
P{s,u) 



where 



Jl = „ „ v . / ^ ^ r Ux, (A.6) 

7(1 - e~^ N ) J 1/N x(l-x) y ' 

1 [ 1/N (e 2yNx - 1)(1 - e-^ Nx ) , 

■h = —. / 7T — \ ~ dx ■ AJ 

7(1 - e t ) J Q x(l -x) 
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A. 2 Branching process theory 

According to Barton (1995), the probability that a beneficial mutation reaches 
fixation in a genotype with genetic background i follows from iterating the 
following set of equations: 

oo 

(1-P i;t _ 1 ) = ^^(l-^, (A.8) 

1=0 

where Pi^ is the probability of fixation of an allele that is present in a single 
copy in site i in generation t, Wi.j denotes the probability that an allele in 
site i contributes with j offsprings to the next generation, and 

P* t = J2Mi, k P k ,t (A.9) 

k 

is the probability that an allele in background i at time t — 1 would be fixed, 
given that at time t it is passed to one offspring. The quantity M^ k gives 
the probability that an offspring from a parent at background i will be at 
background k. If the distribution of offspring is given by a Poisson distribution 
with mean (1 + Sj), i.e., 

W^ = ^-±^-e^ 1+ '*\ (A.10) 

then Eq. (A.8) becomes 

(1 - P iit _!) = exp [-(1 + Si )P* t ] . (A.11) 

The fixation probabilities correspond to the solution of Eq. (A.8) obtained in 
the limit t — > oo. 

In our model, we consider only two distinct classes of genotypes, and mu- 
tations occur only from advantageous mutant to wild type. Because the wild 
type has a mean number of offspring per generation exactly equal to 1, the 
probability of fixation of wild type sequences in the branching-process formu- 
lation vanishes (Harris, 1963). Therefore, the probability of fixation P(s,u) is 
given by the solution of the equation: 

P(s, u) = 1 - exp [-(1 + s)(l - u)P{s, u)} . (A.12) 

Unfortunately, it is not possible to derive a closed-form solution to this equa- 
tion. P(s,u) has to be determined numerically from iterating Eq. (A.12). 

Since the relevant parameter seems to be 7 = (1 + s)(l — u) — 1, the above 
equation can be written as 



P(s, u) = 1 - exp [-(7 + l)P(s, u)} . (A.13) 
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1!) 



Although a derivation can not be performed when n beneficial mutants are 
considered, we found empirically that in this case 



P(s, u) = 1 — exp 



+ 1) 



P{s,u) 



(A.14) 



This expression is used to compare numerical results to the branching process 
theory prediction in Fig. 5. 



